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Abstract 

We present a computation of one- and two-jet cross sections in polarized photon- 
hadron collisions, which is accurate to next-to-leading order in QCD. Our results can be 
used to compute photoproduction cross sections in electron-proton scattering. To this 
purpose, we investigate the structure of the polarized Weizsacker-Williams function, 
where we include a universal, non-logarithmic term, neglected in the literature. We 
construct a Monte Carlo code, within the framework of the subtraction method, and 
we use it to study the phenomenology of jet production in the energy range relevant 
to HERA. In particular, we investigate the perturbative stability of our results, and we 
discuss the possibility of constraining polarized parton densities of the proton and the 
photon using jet data. 
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The last decade has seen an important advance in our understanding of polarized nucleon 
structure functions as a result of the analysis of deep inelastic scattering (DIS) data |]T|. Unfor- 
tunately, the use of DIS data alone does not allow an accurate determination of the polarized 
parton densities. This is true in particular for the gluon, since this quantity contributes to 
DIS in leading order (LO) only via the Q^-dependence of the spin asymmetry {A'^), which 
could not be thoroughly studied experimentally so far. In the near future, the improved ca- 
pabilities of handling polarized beams will allow high-energy colliders to operate in polarized 
mode, thus giving the possibility to study polarized scattering in a way complementary to 
DIS. At the RHIC collider at BNL, the first proton-proton polarized collisions are expected 
in two years from now. The possibility that the electron-proton HERA collider at DESY will 
run with polarized beams has also been considered for a long time [0 . 

At variance with DIS, collider physics offers a relatively large number of processes whose 
dependence upon the gluon density is dominant already at LO. The study of these processes 
is therefore crucial in order to measure this density in a direct way. Among the various high 
transverse momentum reactions with a sizeable gluonic contribution, jet production is an 
obvious candidate, because of the large rates. At HERA, the largest cross sections will be 
obtained by retaining those events where the electron is scattered at very small angles with 
respect to the beam line. In this way, the photon exchanged between the electron and the 
proton is almost on-shell, and the underlying dynamics is that of a photoproduction process. 

As was shown in ref. 0], a polarized version of the HERA collider with ~ 300 GeV 
would be a very promising and useful facility for studying polarized photoproduction reactions. 
In particular, jet production cross sections display a strong sensitivity to the polarized gluon 
distribution of the proton. However, the analysis of ref. is based on a lowest order QCD 
computation. From unpolarized physics it is well known that, in order to describe the data 
in a satisfactory way, the calculation of next-to-leading order (NLO) QCD corrections to the 
jet cross sections is mandatory. 

The computation of the NLO QCD corrections to one- and two-jet production in polarized 
photon-hadron collisions is the purpose of this paper. As is well known, in QCD photoproduc- 
tion cross sections are written as the sum of two terms (point-like and hadronic components), 
neither of which is physically meaningful by itself. The hadronic part, in which the photon 
behaves like a hadron, can be accounted for by using the available results for jet production 
in polarized hadron-hadron scattering [Q. The NLO QCD corrections for the point-like part 
are presented in this paper for the first time. Such a calculation needs the one-loop 2 — > 2 and 
tree-level 2 — > 3 polarized amplitudes as input (one of the particles in the initial state being 
the photon and the other a quark or gluon). These amplitudes are already known ^j. The 
next step is to implement them in a Monte Carlo code which allows the calculation of any 
three-parton infrared-safe observable. In order to do this, we will use the general algorithm, 
based on the subtraction method, presented in ref. [^]. Following the same strategy as already 
adopted in the case of polarized hadron-hadron collisions 0], we build a modified version of 
the Monte Carlo code of ref. ||^, which deals with unpolarized photon-hadron collisions. The 
reader can find further details in ref. 

In photoproduction processes initiated by an electron, the electron can be considered to 
be equivalent to a beam of real photons, whose distribution in energy (Weizsacker-Williams 
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function 0) can be computed theoretically. The polarized electron-hadron cross section reads 



Vmax 



dAaMKe,KH) = / dyAf\'^{y)dAa^H{yKe,KH), (1) 

where Kg and Kh are the momenta of the incoming electron and hadron respectively, and 
< ymin < Umax < 1 are fixcd by kinematical boundary conditions or experimental cuts. The 
quantities Aa appearing in eq. (|I]) are defined in terms of cross sections cr(Ai, A2) for incoming 
particles of definite helicities 

dAa = ^[d(x{+, +) + d(x{-, -) - da{+, -) - da{-, +)). (2) 

The unpolarized electron-hadron cross section can be obtained from eq. (|l]) with the formal 
substitutions 

Aa-^a, A/('=) /(^). (3) 

As mentioned before, the polarized and unpolarized Weizsacker-Williams functions (A/^*^) 
and f!f \ respectively) can be computed. As far as the unpolarized case is concerned, we use 
the form M 



(.(e)/ ^ Q'fi 
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(4) 



Notice that eq. (^) contains a non-logarithmic term with a singular behaviour for y ^ 0. In 
ref. [10 1 this term was shown to give non- negligible contributions (of about 7% and 5% for 
= 0.01 GeV^ and = 4 GeV^, respectively) in the case of unpolarized jet production 
at HERA. One expects an analogous non-logarithmic term to appear also in the polarized 
Weizsacker-Williams function. To the best of our knowledge, however, such a term has never 
been presented in the literature. We therefore computed A/^*^^ from scratch. We closely 
follow the procedure of ref. [^, in which the interested reader will find further details. We 
consider the generic production process 

e{p; n) + a{k; s) — > e{p') + X, q = p-p', (5) 

where a is a parton whose momentum k, which is kept off-shell for the time being, will 
eventually be put on-shell, and n, s are spin vectors. In the case of polarized scattering, 
we only need the antisymmetric part of the partonic and leptonic tensors (see for example 
ref. Q): 



5(i(g^ k ■ q)s„ + g2{q^, k ■ q) (s„- ^^k^ 



(6) 



Tr = ^ = -2^mee^^^'^npg., (7) 

where gi and g2 are the usual structure functions relevant to polarized collisions. The spin 
vectors for the incoming electron and parton can be written as 
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where Nn and Ng are normalization factors such that \n^\ = = 1. The spin vectors satisfy 
the transversahty conditions n ■ p = and s ■ k = 0. The electron-parton cross section follows 
from direct computation]^ {y = {k ■ q)/{k ■ p), k"^ — > 0): 

1 e^W^'TA,. d'p' 



dAa^aip, k) 



Ak-p {27r)^2E' 

1 - (1 ^ 2m2?/2 



an 
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—n dq dy. 

Ay k ■ p 



(9) 

yq. - - 

Notice that the terms proportional to g2 drop from this equation. This expression can 
be related to the polarized cross section for real photon-parton scattering, which can be 
computed by convoluting the parton tensor of eq. (|) with the antisymmetric part of the 
photon polarization density matrix, 



2 ' ■ 2^|g2| 

where is the photon polarization vector and is its spin vector 

t^. = Nt [q, 



We get (g^ ^ 0) 



Combining eqs. 
ref. [1), we get 



Acr^fa{q,k) 



1 



^'^^"''^^"4A;-g"^ ""^^ 4k ■ q 
and integrating over in the accessible kinematical range (see 



5-1(0, k ■ q) 



(10) 
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and 



dAaea{p, k) = Aa^^iyp, k) A&\y) dy 



where 



27r 
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1)! ,0, (-1-1^ 



(13) 



(14) 



We therefore find that also in the case of polarized scattering a non-logarithmic term is present. 
This term, at variance with the unpolarized case, is not singular for y ^ 0. However, exactly 
as in the unpolarized case, its behaviour for small y is equal to the behaviour of the term that 
multiplies the logarithm. 

We can now go back to eq. (|l|), and compute dAa^n- Using the factorization theorems of 
QCD, we can write 

dAa^niK^, Kh) = dAa^^t\K^, Kh) + dAa';t{K^, Kh), (15) 

where 

dAa^H''\Ky, Kh) = Y1 J dxAf^"\x, iJ,'„)dAa^j{K^, xKh, asifi'j^), fi'j,, fi'^, fJ-y), (16) 

j 

dAa^^^^iK^Kn) = J dxdyAf^'\x, fi,)Af^''\y, f,';,) 

X dAaij{xK.,, yKn, as{^'^), fi'^, /x'^, fi^) . (17) 



^We keep only the first term in the expansion of gi and 52 in series of q^: the remaining terms give 
non-factorizable corrections which are suppressed by powers of Q^/M^. 
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We remind the reader that the two terms in the RHS of eq. (|T5]) (denoted as point-hke and 
hadronic components) are not separately well defined in perturbative QCD beyond leading 
order: only their sum is physically meaningful. The polarized parton distribution functions 
A//^'' and Af^"'^ are defined, as usual, in terms of densities for partons of definite helicity in 
hadrons of definite helicity {A = -ff, 7): 

The subtracted partonic cross sections dAcr^j and dA&ij are finite at any order in QCD and, 
for the purpose of this paper, have been evaluated to NLO (aem«| and respectively). 
Notice that while the dependence upon the scales /i'^,/^^ and cancels (up to NNLO 

terms) in the point-like and hadronic components separately, the dependence upon /x^ in 
the point-like component is compensated by a corresponding dependence in the hadronic 
component. We finally point out that equations completely similar to eqs. (|T3p-(|I^) also hold 
for unpolarized scattering. 

In what follows, we will present phenomenological predictions relevant to HERA physics. 
To this purpose, we will evaluate both the polarized and unpolarized cross sections. Since 
one of the major goals of high-energy colliders in the polarized mode will be that of studying 
the polarized parton densities, we will make a definite choice for the unpolarized densities 
throughout the paper, and study the dependence of our predictions for the asymmetries upon 
the polarized densities. Namely, we will adopt the set MRSA' []T2[ for the unpolarized proton. 



and the set GRV-HO [|T3l for the unpolarized photon. The value of Aqcd associated to MRSA' 
(A^^ = 152 MeV) is well below the current world average; however, it is closer to the value 
of lyqcD associated to GRV-HO and to all the polarized sets we are going to use. The use of 
more modern proton sets would only affect the absolute normalization of the cross sections, 
without any sizeable change in the shape of the distributions. 

As far as the polarized densities in the proton are concerned, DIS data available at present 
leave the gluon density largely unconstrained, while they determine the quark densities to a 
reasonable extent. Several parametrizations are therefore available, whose gluon densities are 
very different from each other. We will consider the following sets: GRSV STD and GRSV 
MAXG 0, DSSl, DSS2 and DSS3 p[, and GS-C Il6| (see ref. for a discussion about 



this choice and a comparison between the various sets). GRSV STD will be the default set 
when studying the perturbative stability of our results. 

The polarized densities in the photon are completely unmeasured so far, and models for 
them have to be invoked. To obtain a realistic estimate for the theoretical uncertainties due 



to these densities, we use the two very different scenarios considered in ref. [jl^, assuming 
"maximal" (A/(T)(a;, /if) = /(^)(x,/if)) or "minimal" {Af^^\x, ii]) = 0) saturation of the 
positivity constraint^ |A/*^'^)| < f^"'^ at the input scale /i/ (C(0.5 GeV)) for the QCD evolution 
(here, the unpolarized densities are those of the set GRV-HO). The corresponding sets are 
denoted as SV MAX 7 and SV MIN 7 respectively. The former will be our default choice 
when studying perturbative stability. 

In order to compute the NLO results consistently, we will use the NLO-evolved version 
of the previously mentioned parton distribution^. We point out that these NLO-evolved 

^Strictly speaking, this bound is violated in QCD, and the deviations may become sizeable at small scales. 
For a discussion on this point, see for example ref. jl^ 

•^In fact, for most of the sets considered here, there exists a LO-evolved version as well. 
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sets will also be used to compute the Bom results; indeed, as observed in ref. 0], the use of 
LO-evolved parton distributions could introduce differences between the LO and NLO cross 
sections that are much larger than those that can be expected on the basis of perturbative 
considerations. 

Since this paper presents the first NLO calculation of jet cross sections in polarized photon- 
hadron collisions, it is mandatory to assess the effect of the radiative corrections. Lacking 
an NNLO computation, the most reliable estimate of the theoretical uncertainty affecting 
our results comes from the study of the dependence upon the mass scales entering eqs. ( [16|) 
and (|13), at fixed values of the other input parameters. Unless otherwise specified, we will 
set all the scales equal to a common value /i, and vary fi in the range < fj, < 2/io. Here 
Ho is the default scale, which we set equal to half of the total transverse energy of the event; 
other sensible choices for fiQ give completely similar results. From eq. (||), it is clear that all 
the dynamical information on the ep process in the Weizsacker- Williams approximation are 
contained in the jp cross section. Therefore, we will study the scale dependence of our results 
using monochromatic photon-proton collisions at E^^'> = 268 GeV; the scale dependence at 
lower centre-of-mass energies, and in electron-proton collisions, is comparable or milder. 

We start by considering the single-inclusive jet pseudorapidity, requiring pj, > 10 GeV, 
where Pt is the transverse momentum of the jet. The jet has been defined using the cone 
algorithm with R = 1. The results for polarized and unpolarized cross sections are presented 
in fig. |1]. The size of the radiative corrections is larger in the case of unpolarized collisions; 
in the case of polarized collisions, the contributions from the various partonic subprocesses 
do not have the same sign, and cancellations occur. In both cases, the variation of the scale 
induces a variation of the cross section of the order of 10% over most of the t] range considered. 
The scale dependence is strongly reduced when going from Born (boxes and circles) to NLO 
(dashed and dotted histograms) results. However, this is not true in the region of negative 
?7's, where the scale dependence of the Born cross section is very small, and smaller than that 
of the NLO results. This could be the signal of a failure of the perturbative expansion, which 
would be extremely relevant from the phenomenological point of view, since it is precisely 
from this region that one hopes to extract information on the gluon density in the proton . 
In order to understand this point, we have studied the scale dependence of our results more 
carefully, by varying separately all the scales that enter eqs. (|16D and (0). At the Born level, 
the hadronic component has a relatively mild /i'^ dependence, a sizeable fi^ dependence, and a 
very large /i'^ dependence. Although the /i^ and /i'^ variations tend to compensate each other, 
the latter is dominant, and the overall result is a large scale dependence, which is precisely 
the effect we see in the positive rj region in fig. |^, where the hadronic component is dominant. 
On the other hand, the point-like component displays a dependence upon /i'^ and fi'^ of about 
the same size, but opposite in sign (at the Born level, the point-like component does not 
depend upon /i^), and a simultaneous variation of the two scales leads to an almost perfect 
cancellation of the effects induced by them. This is the origin of the small scale dependence 
of the Born result that we see in the negative 1] region in fig. |T]. When going to NLO, in the 
hadronic contribution the /i'^ dependence basically reduces to zero, and the ^'^ dependence 
is sizeably reduced; this is what we expect, as remarked after eq. ([181) . The fj,^ dependence is 
also reduced, although it remains sizeable; indeed, the effects of the variation of fi-y are only 
partially cancelled in the hadronic component, since the complete cancellation (up to NNLO 
terms) takes place when the point-like component is included. In fact, we can observe this 
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Figure 1: Scale dependence of the single-inclusive jet pseudorapidity in polarized 
(left) and unpolarized (right) collisions. 



cancellation by looking at the /x^ dependence of the point-like component at NLO, whose 
size is almost the same as that of the hadronic component, but opposite in sign. Finally, the 
^'jj and /i'^ dependences of the point-like component are largely reduced with respect to the 
corresponding dependences at LO. In conclusion, the scale dependence of the rj distribution is 
consistent with what we expect from perturbative considerations. The small scale dependence 
of the Born result in the negative r] region is only due to an incidental cancellation, which 
takes place when fixing all the scales to the same value. A more sophisticated treatment of 
the scale dependence would result in a scale dependence larger at the Born level than at the 
NLO level. 

We also studied several double differential observables, by considering the two leading jets 
of each event. As an example, we present in fig. H the invariant mass of the pair. We require 
the hardest (next-to-hardest) jet to have transverse momentum larger than 15 GeV (10 GeV). 
We do not impose any t] cuts, in order to maximize the scale dependence (as shown in fig. |l], 
the scale dependence is larger for the smallest and largest accessible rj values, since they 
correspond to small transverse momenta). Therefore, in the case of rj cuts simulating a realistic 
geometrical acceptance, the scale dependence of our results would be smaller. We see that 
the case of the invariant mass is pretty similar to the case of single-inclusive pseudorapidity. 
The size of radiative corrections is larger in the case of unpolarized scattering, and the scale 
dependence of the polarized result is comparable to that of the unpolarized one. A reduction 
of the scale dependence can be observed when going from LO to NLO, except for the region 
close to the threshold, where there is no contribution at the Born level (due to the asymmetric 
transverse- momentum cuts), and the NLO result is effectively a LO one. 

For all the single-inclusive and double-differential jet observables that we have studied, 
the same pattern as outlined above is reproduced. The size of radiative corrections is usually 
smaller in the polarized case than in the unpolarized case. The inclusion of NLO terms reduces 
the magnitude of the scale dependence with respect to the LO results, in both the polarized 
and unpolarized cases. When considering double-differential cross sections, special care has 
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Figure 2: Scale dependence of the invariant mass of the jet- jet pair in polarized 
(left) and unpolarized (right) collisions. 



to be taken in the case where the cuts on the transverse momenta of the two leading jets 
are set to the same value, since there are regions in the phase space where the perturbative 
computation fails to give a sensible prediction. This situation has been discussed at length 
in ref. |TU1 for unpolarized jet production. We did not find any major difference in the case 



of polarized scattering. We also investigated the effect of changing the jet definition, by 
considering the kj. algorithm proposed in ref. [|r^ . This basically only amounts to a change in 
the normalization, while the effects due to the variation of the scales are of the same size as 
those studied above. Finally, the perturbative stability of the results obtained with different 
parton densities is identical to the one presented here. 

We now turn to the problem of studying the dependence of our results upon the available 
proton and photon polarized parton densities. To this end, we will consider electron-proton 
collisions in the Weizsacker-Williams approximation, setting E^^'^ = 300 GeV, ymin = 0.2, 
Vmax = 0.8, and = 0.01 GeV^ (see eqs. (|I]), (D and (|I3|)). All the scales will be set to the 
default value. In fig. ^ we present the results for the asymmetry 

daep/dr) 

where r] is the pseudorapidity of the single-inclusive jet. A cut of pr > 10 GeV has been 
applied. On the LHS of the figure, we have chosen GRSV STD as the polarized proton set, 
and we evaluated the results for both the polarized photon sets we consider in this paper. The 
Born and NLO results are both shown. We can see that in the large (positive) rj region the 
difference induced by the choice of the two photon sets is extremely large. On the other hand, 
towards negative 7] values this difference tends to vanish. This is because in that region the 
point-like component, which does not depend upon photon densities, is the dominant one. We 
can also observe that in the positive t] region there is a very small difference between the NLO 
and LO results, while for negative 77's the radiative corrections are positive and reduce the 
asymmetry considerably. It is worth noticing that, as can be seen from fig. this reduction 
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Figure 3: Asymmetries versus pseudorapidity in single-inclusive jet production, for 
electron-proton collisions at Ecm = 300 GeV. The predictions for several photon (left) 
and proton (right) densities are shown. 



is mainly due to the increase of the unpolarized cross section at NLO with respect to LO, 
while the polarized result is basically unchanged. On the RHS of fig. we show the curves 
obtained by fixing the polarized photon set to SV MAX 7, and by considering the various 
polarized proton sets. As expected, the largest differences can be seen at negative rj values, 
where theoretical predictions can vary for about one order of magnitude. The shape of the 
asymmetries is quite similar for all the sets, except for GS-C; this reflects the fact that the 
shape of the GS-C gluon is rather different from the shape of the gluons of the other sets. 
The same effect has indeed been observed in polarized pp collisions @] . 

On the LHS of fig. ^ we also show (dotted lines) an estimate of the minimum value of the 
asymmetry observable at HERA. This is calculated using the following formula 

(A) ■ = — — , ^ (20) 

where C is the integrated luminosity, Pe(p) is the polarization of the electron (proton) beam, 
and the factor e < 1 accounts for experimental efficiencies; is the unpolarized cross section 
integrated over a range in rapidity (A77), which has been chosen, for the present case, equal to 
0.2. In fig. ^ we used C = 100 pb~^, Pe = Pp = 1 and e = 1. It follows that, if high luminosity 
will be collected, it will be possible to get information on the polarized parton densities in 
the proton, given the fact that realistic values for the polarization of the beams will be such 
that PePp — 0.5. Notice that the cross section corresponding to sets such as DSS3 or GS-C 
will be hardly measurable at HERA, regardless of the luminosity collected. Of course, the 
bin size At] can be enlarged; this would result in a smaller observable asymmetry, at the 
price of a loss in resolution. As far as the polarized photon densities are concerned, if the 
"real" densities are similar to those of the set SV MIN 7, it will be extremely hard to even 
get the experimental evidence of a hadronic contribution to the polarized cross section. On 
the other hand, a set like SV MAX 7 appears to give measurable cross sections, but this 
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conclusion strongly depends upon the polarized proton densities. As observed in ref. 0, 
sensible information on the polarized photon densities will be obtained only after a more 
precise knowledge of the proton densities will be availabe. At this point, we have to comment 
upon the effect of the non-logarithmic term in the polarized Weizsacker- Williams function. It 
turns out that, for the single- inclusive distribution we are studying here, this term reduces the 
absolute value of the cross section by a factor of about 6%. Since the non-logarithmic term 
in the unpolarized case also reduces the cross section by a factor of about 7% [0], it follows 
that the asymmetry is to a good extent independent from the presence of non-logarithmic 
terms in the Weizsacker- Williams functions. Finally, if a different jet definition is adopted, 
the asymmetries change for less than 1%; this feature has already been observed in the case 
of pp collisions ^ . 

In principle, it is conceivable to study asymmetries also for monochromatic photon-proton 
collisions. From the experimental point of view, this would correspond to selecting events with 
the Weizsacker-Williams y (that is, the fraction of electron energy carried away by the pho- 
ton) in a given narrow range (which corresponds to taking y^m and Umax very close to each 
other in eq. (0)). If we compute the asymmetry of eq. (|19[) with {A)a^p instead of (A)(Jep, the 
values we obtain are larger than those shown in fig. |[ For Ej^^^ = 134 GeV, the asymmetry 
can be twice as large as that obtained using eq. (jT^). Unfortunately, this is not sufficient. 
Indeed, what is actually measured at HERA is an electron-proton cross section. Therefore, 
to get a physically observable quantity one should multiply the 7p asymmetry by the ratio 
of the integral of the polarized Weizsacker-Williams function in the range i/min < 1/ < Umax 
over the same integral calculated for the unpolarized function. This ratio tends to zero for 
y mini Umax ^ {ymax — Umin IS kept fixcd), thus disfavouriug the measurement of the asym- 
metry at small photon-proton centre-of-mass energies. 

In conclusion, we reported the first calculation of jet cross sections in polarized photon- 
hadron collisions, which is accurate to NLO in perturbative QCD. For all the observables 
considered, it has been found that the size of the radiative corrections is moderate, and that 
the scale dependence is smaller than that of the LO result. The inclusion of the NLO terms 
reduces the size of the asymmetries in the pseudorapidity region where the contamination 
from the hadronic photon contribution is minimal. In order to get information on the po- 
larized gluon densities in the proton and in the photon, a very high luminosity must be 
collected at HERA. We also calculated the universal, non-logarithmic term in the polarized 
Weizsacker-Williams function. Its contribution to the cross sections is small with respect to 
the logarithmic term, analogously to what happens in the unpolarized case. 
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